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We consider slow, steady transport for the normal state of the superconductor La2Cu04+s in 
a one-dimensional geometry, with surface fluxes sufficiently general to permit oxygen to be driven 
into the sample ( "loaded" ) either by electrochemical means or by high oxygen partial pressure. We 
include the bulk reaction O— +0 2 ~ + 2h, where neutral atoms (a) go into ions (i) and holes (h). 
This system is a mixed ionic electronic conductor (MIEC). For slow, steady transport, the transport 
equations simplify because the bulk reaction rate density r and the bulk loading rates dtn then are 
, uniform in space and time. All three fluxes j must be specified at each surface, which for a uniform 

current density J corresponds to five independent fluxes. These fluxes generate two types of static 
modes at each surface and a bulk response with a voltage profile that varies quadratically in space, 
£N) , characterized by J and the total oxygen flux jo (neutral plus ion) at each surface. One type of 

surface mode is associated with electrical screening; the other type is associated both with diffusion 
and drift, and with chemical reaction (the diffusion-reaction mode). The diffusion-reaction mode is 
accompanied by changes in the chemical potentials jj,, and by reactions and fluxes, but it neither 
carries current (J = 0) nor loads the system chemically (jo = 0). Generation of the diffusion- 
reaction mode may explain the phenomenon of "turbulence in the voltage" often observed near 
the electrodes of MIEC's. Within the bulk, the local fluxes satisfy a relation that is independent 
of the applied fluxes. As a consequence, the bulk response alone cannot match arbitrary values 
for the five independent input fluxes; matching occurs by generating appropriate amounts of the 
diffusion-reaction mode at each surface. The bulk response is completely responsible for steady-state 
loading and typically possesses a voltage profile that varies quadratically in space, as for the lead- 
acid cell. Seven macroscopic parameters (three d/j,/dn's, three diffusion constants, and a reaction 
rate constant) characterize the theory. 
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PACS numbers: 74.72.-h, 82.33.Pt, 82.45.Xy 



I. INTRODUCTION 



Mixed ionic electronic conductors (MIECs) are often "loaded" with a specific atom, in either neutral or ionic form. 
In La2Cu04+a, atomic or ionic oxygen is loaded into a sample in order to produce a concentration that is more 
favorable to superconductivity^^ An idealization of such an experiment would consider a rod of the material in a one- 
q \ dimensional geometry where fluxes of neutral atoms a (O), ions i (O 2- ), and holes h slowly and steadily enter or leave 
each end (at x = and x = C). Thus, in addition to the material parameters of the system (three thermodynamic, 
three transport, and one reaction), steady state transport requires specification of six fluxes, subject to the condition 
that the same net current density J enters and leaves, making a total of five independent fluxes at the two surfaces. 
Hence, although the rod's net oxygen content (neutral plus ion) increases with time, its net electric charge remains 
constant. Once the neutrals or ions enter the rod, they may recombine via the reaction O^O 2- + 2h. 



A. Chemical Loading and Chemical Reactions 



Slow steady transport with loading but without chemical reactions has already been studied for the multiple charge- 
carrier systems of lead-acid cell^ and La2Cu04+5-- Because the diffusivities of the charge-carriers differ, the voltage 
profile within the bulk has a component that varies quadratically in space, the details depending on the fluxes at the 
surface. These systems support surface modes of only one type, which correspond to electrical screening. 

Within the context of semiconductors, slow steady transport with chemical reactions (electron-hole recombination) 
but without loading has also been studied^ Because of recombination (a form of chemical reaction) there is an additional 
surface mode, but no quadratic variation in space of the voltage profile. The additional mode is associated with 
diffusion and drift, as well as with recombination, and varies exponentially in space. Although its amplitude is 
dependent on these fluxes, its structure is independent of the fluxes at the surface. 

The present work considers slow steady transport with both loading and chemical reactions. The example again is 
La2Cu04+a, since whatever the excess oxygen content 5, there is a way to force oxygen to enter such that reactions 
must take place to cause equilibration. For example, if all the oxygen enters the bulk as ions (which can be detected 
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sample of I^CuO^^ 



FIG. 1: Rod of La2Cu04+5 with two reference electrodes attached near its ends. 



by the conductivity of the sample 1 ), then by employing a high pressure atmosphere of molecular oxygen one can 
expect atomic oxygen to enter the sample, and then convert by reaction to ionic oxygen. Typically, when oxygen is 
added the fraction that goes into atomic bulk states is neither zero nor unity, so that bulk reactions are needed to 
cause bulk equilibration.^ Therefore, we consider an experiment where the surface fluxes of atomic and ionic oxygen 
are completely arbitrary. 

In practice, the specific values of the surface fluxes may or may not be known in a given experiment. Assuming 
that they are known, in this work we show how one may obtain the response of the MIEC, thus giving the rate of 
loading of O and O 2- , the potential, the densities, and the fluxes throughout the sample. Our analysis is valid in 
the limit of low fluxes, where the equations can be linearized, and for slow steady fluxes. In principle, the fluxes can 
vary slowly in time. (Slow is loosely defined relative to time it takes for the component with the smallest diffusivity 
to diffuse across the sample.) Corresponding to the six fluxes specified at the boundaries there must be: at a given 
surface, a combination of fluxes that specify the amplitude of the exponentially-varying diffusion-reaction mode at 
that surface; for the bulk, four independent flux combinations that specify the bulk response, which gives a uniform 
chemical loading within the sample and a uniform current flow through the sample. To our knowledge, this is the 
first work on any system to consider diffusion and drift, bulk reactions, and surface fluxes general enough to include 
chemical loading. 



B. "Turbulence" at the Interface 



Voltage measurement in the area of MIECS is not trivial. The review by Kudo and Fueki^ notes that "In actual 
measurements, however, there is often turbulence of the potential distribution in the vicinity of the electrodes with 
which the electric field is applied." The authors do not define "turbulence" in any operational sense; it may refer to a 
complex variation in space or in time. It has become conventional for experiments in this area to place the reference 
electrodes away from the contact between the MIEC and the adjacent conductors (electronic or ionic), presumably 
to avoid such "turbulence." See Figure 1. 

If turbulence refers only to a complex variation in space, then it could be due to surface modes generated near 
the interfaced There are at least two such modes in systems of this sort: the usual Debye-like screening mode 
and the above-described diffusion-reaction mode (in more complex systems, there can be multiple diffusion-reaction 
modes). The reference electrodes measure "voltages" re / proportional to electrochemical potentials fx (4>h = fih/e 
for hole conduction, and <f>i = — jli/2e for divalent anion conduction). Corresponding to the "voltages" are the "fields" 
Eh.i = —d x 4>h,i- Since the screening mode involves no change in the electrochemical potentials, the reference electrodes 
do not respond to the screening mode. On the other hand, the electrodes do respond to the diffusion-reaction mode. 



C. Outline of Paper 



Sect. II employs the methods of irreversible thermodynamics, including the effects of chemical reactions, to obtain 
the transport equations for mixed ionics, including the effects of reactions. Sect. Ill obtains the steady-state surface 
modes. Sect. IV obtains the steady-state bulk response. Sect.V discusses how to extract the parameters of the system 
from measured quantities. Sect. VI presents our conclusions. 
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II. IRREVERSIBLE THERMODYNAMICS 



Consider, in one-dimension, a uniform sample of I^CuG^+a-, as in Fig.l. We take the carriers to be holes and ions 
(O 2- ) produced by neutral oxygen atoms O via the reaction 

0^0 2 -+2h (1) 

with reaction rate density r. 



A. General Considerations 



Let u denote the energy density, T the temperature, s the entropy density, p^, Pu and p a the hole h, ion i, 
and neutral atom a electrochemical potentials, and n^, Uj, and n a the corresponding number densities. Then the 
fundamental thermodynamic differential for this system is 

du = Tds + jlhdnh + fiidm + p a dn a . (2) 

With /i/j, fii, and p a the chemical potentials, and </> the electrical potential, we have 

Aft = f*h + e 

Here 4> satisfies Poisson's equation 

V 2 = 

where e is the dielectric constant, a multiple of the permittivity of free space Eq, and the charge density 

p = e(-2m + n h ). (5) 
The conservation laws for this system are, in one-dimension (x) 

d t u + 8 x f = 0, d t s + d x f = ^ > 0, (6) 



b, p,i = pi- 2e</>, p a = p a . 



p e . . 
- = --{-2m + n h ), 



(3) 
(4) 



d t n h + d x j h = 2r, d t n t + d x f = r, d t n a + d x f = -r. (7) 

Here j u is the energy flux density, j B is the entropy flux density, V is the rate density of heat production (V/T is 
rate density of entropy production), j h is the hole number flux density, j l is the ion number flux density, and j a is 
the neutral atom number flux density. In (JJJ, the signs of r are taken such that if O— >0 2 ~ + 2h, then nh/2 and rii 
increase at the same rate that n a decreases. The fluxes, r, and V are to be determined. 
The time-derivatives of (0 and © (0 lead to 

< V = -dij u + Td x f - p h (2r - d x] h ) - pi(r - d x f) - p a (-r - d x f) 
= -d x (j u - Tf - p h] h - p l3 l - p a f) 

-fd x T - j h d x p h - fd x p,i - j a d x p a - r{2p h + pi - p a ). (8) 

Expressing V as a non-negative quadratic form uniquely requires that: (a) 

r = Tf + p h j h + iMf + p a j a ; (9) 

that (b) 

f = -j,d x T - a sh d x p h - a si d x pi - a sa d x p a , (10) 

where k > is the thermal conductivity; that (c) 

j h = -a hh d x p h - a hi d x pi - a ha d x p a - a hs d, t T, 
f = -a lh d x p h - a ti d x pi - a ia d x p a - a is d{T, 
j a = -a ah d x p h - a ai d x pi - a aa d x p a - a as d l T, (11) 
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where a hh > 0, an > 0, and a aa > 0, (n/T)a h h > o? hs , {n/T)a vi > af s , (n/T)a a a > a 2 as , and a^aa > o? hi , etc. (by 
the Onsager symmetry principle 8 , — onh, a s h — oth s , and a s i — oti S ); and that (d) 

r = -A(2/}/j + p,i - /j, a ) = -\{2u h + Hi - [i a ). (12) 

where the reaction parameter A > is related to the hole and ion lifetimes th and n. In equilibrium r = 0, so l|12|) 
implies that the change in Gibbs free energy be zero, or = 2jj^' + — fJ. a °\ as expected. For slow steady processes, 
r need not be zero but, as we will see, r will take a value consistent with the surface fluxes, and r will be uniform in 
space and in time. 

B. Specific Considerations 

From here on, we neglect any coupling to entropy or temperature, and any off-diagonal coupling. Thus = Oj s = 
ahs = a as = a hi = a ha = a la . Moreover, we employ a h = a hh , a a = a aai and a, t = a u . 
We write the densities in the form 

n = n° + An, (13) 

where n° is the equilibrium value and An is the deviation from equilibrium (Sn will be reserved for terms that vary 
in space, either in the surface solution or in part of the bulk response) . We now linearize r as 

r = -X(2u h + m - Mq ) = _A(2^An, + ^An t - §^An„). (14) 

dn h dni dn a 

Defining the reaction rates w (with dimensions of inverse time) as 

Wh = — , Wi = A—, w a = A- — , (15) 
dn h dni on a 

(fH|) becomes 

r = -(2w h An h + WiAn t - w a An a ). (16) 

Because the w's vary as the product of a thermodynamic derivative and the reaction parameter, they can be thought 
of as "thermo-reaction" parameters. 

Linearizing the fluxes and neglecting the off-diagonal terms, (|11|) becomes 

j h = ~a h d x jl h = -a h d x (fi h + efi) = ~a h (^-d x n h + ed x (f>), (17) 

dnh 

f = -otid x p,i = -aid x (ui - 2ecp) = -Oi(-^d x ni - 2ed x <j)), (18) 

j a = -a a d x fi a = -&ad x u a = -a a ^-^-d x n a . (19) 

dn a 

The charge-carrier conductivities a are related to the a's via 

a h = e 2 a h , er; = Ae 2 ^. (20) 
In terms of the effective electric fields for the charge carriers, 

E h = -d x (f>h, Ei = -d x ct>i, (21) 

we have 

j h = a h E h , j^-la^. (22) 

The diffusivities D are given by 

E>h = a h -^ — , E>i = Oii- — , E) a = a a - — . (23) 
dnh dn t dn a 
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Because the D's vary as the product of a thermodynamic derivative and a transport parameter, they can be thought 
of as "thermo-transport" parameters. Note that the ratio of a w to a D depends only upon reaction and transport, 
the thermodynamic derivatives cancelling. Such ratios occur for the diffusion-reaction surface mode. 
For completeness, we observe that the electric current density is given by 

J = e(j h - 2j% (24) 

and the net oxygen mass flux is given by 

jo=f+j a . (25) 



III. STEADY-STATE SURFACE MODES 



For surface modes we will employ the notation 5n for the deviation of a number density from equilibrium, and 5(f> 
for the deviation of the potential from equilibrium. We also assume that there is a local relation between n and fj,, so 
that © can be employed to relate 5n to 5ft, and 5(f>. Then the linearized form of (0J becomes 

V 2 ^ = -(25m - 5n h ) = -(p^ + 4^)50 + -(2^ft - P±6fi h ). (26) 
e e d[i h oui e dui o\i h 



A. Screening Surface Mode 

A steady-state (dn/dt = 0) solution of the equations occurs for no shift in the electrochemical potentials (5 ft = 0). 
In this case local equilibrium holds, so there are no fluxes and the recombination rate r = 0. As noted above, this 
mode cannot be measured by electrodes sensitive to individual electrochemical potentials, 9 although the associated 
density changes might be susceptible to other types of analysis (optical, chemical, or otherwise). Setting all the 
electrochemical potentials to zero in Q26JI then gives 

V * 5(t >= e -(pL + 4p)5ct>. (27) 



This has solution 



1 e 2 ,dnh . dm 



I 2 e y dfih dm 



Aexp(~x/l)+Bexp(x/l), = — (t^+ 4 7T^)- ( 28 ) 



This screening mode has screening, or Debye, length I given by l|28(l . This mode depends only upon equilibrium 
properties of the system. Note that the neutral atoms do not participate at all. 



B. Diffusion and Reaction Surface Mode 



Another steady-state (dn/dt — 0) solution is obtained by considering the case where r ^ 0. Comparison of the 
three terms in the continuity equations JIJ) gives the conditions 

j h /2 = f = -f. (29) 

Note that H29J) implies J = (no current flow) and jo = j l + j a = (no net oxygen flow). By the flux equations Ulljl. 
the variations in electrochemical potential then satisfy 

-ah5p,h/2 = -aiSpi = a a 5jl a . (30) 

Thus H12fl can be rewritten as 

r « -2\a h 5ji h {— + + — ). (31) 

Using this in the steady-state version of (JJJ for j h then gives 

-cthdlfih = ~4:Xa h 5fL h (— + — + — !— ). (32) 
ah 4ai 4a a 
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This is solved by 



= Aexp(-x/L)+Bexp(x/L), -L = 4A(— + + -^). (33) 

i> a/i 4a, 4a a 



This diffusion-reaction mode has diffusion-reaction length L given by (|33|l . All three mobile species contribute to this 
mode. Since the a's are related to diffusion, this mode depends both upon diffusion (a) and reaction (A). The faster 
the diffusion (i.e. the larger the a), the larger the L; the faster the reaction (i.e. the larger the A), the shorter the 
L. This qualitative dependence on diffusion and reaction is as expected. A similar mode occurs for semiconductors, 
where the chemical reaction involves the recombination of electrons and holespi 

If any of these mobile species has a very slow rate of transport, corresponding to a small a, then L is very short. 
As a consequence, within a short distance of the surface the system can adjust from surface-determined boundary 
conditions to the bulk values. 

For this mode, if 5p,h = (C/a/j) exp(— x/L), then by (|30|l 

C C C 

Sfi h = — exp(-x/L), Sfii = - — exp(-x/L), 8fj, a = -- — exp(-x/L). (34) 
oth 2a t . 2a a 

By (I26|) the potential 5(f> is given by 

/ 1 1 \xj. e 1 dn h 1 

(T2 _ 72^ = — a )Cexp(-x/L). (35) 

This equation relates A of Q to C of 

Since the potential <5(/> of (|33|l is non-zero for the diffusion-reaction mode, this mode also involves the electric field. 

Finally, note that associated with any steady-state surface mode, where d t n — 0, there is no deposition of material 
near the surface. For the present mode, where there are both fluxes and chemical reactions, whatever component is 
produced by chemical reactions (r 7^ 0) is taken up by nonuniform flux (d x j =/= 0). 

If jo = j l + j a is non-zero for the bulk response, the diffusion-reaction mode (with jo non-zero but j % and j a 
individually non-zero) has the important property that, when added to the bulk mode, it can change the flux ratio 
j l + j a ) on moving from the surface to a few L within the bulk. 



IV. STEADY BULK TRANSPORT RESPONSE 



There are two types of bulk transport response, according to whether or not the system is being chemically loaded. 



A. Steady Bulk Transport Response — No Chemical Loading 



One solution of the steady-state equations (where dtnh = 0, etc.) occurs for Sn^ = 5n. L = 5n a = 0, but 8cj> — 
—Ex + D, corresponding to a uniform shift in the electrical potential and a uniform field E x — E. This leads to 
constant partial currents j h = cth,E, j l = —2a.iE, j a = 0, and corresponds to no chemical loading, although there is 
a net oxygen flux that crosses the system. The solution we obtain in the next section is sufficiently general that it 
includes this case, which will serve as a check on the results of the next section. 



B. Steady-state Bulk Flow with Chemical Loading, so d t n 7^ 



Consider a situation with two electrode surfaces, one at x — and the other at x = C. Let there be slow steady 
flow, and at the surfaces let the ion, atom, and hole fluxes j l , j a and j h be specified, subject to equal net electric 
current densities J at the surfaces (so there is no electrical charging). However, there can be unequal net oxygen fluxes 
j l + j a at the surfaces, so a net amount of oxygen can be loaded into the system (mass charging). Thus, we specify 
the fluxes at the boundaries, and within the system we must determine the potential, the fluxes, and the densities. 
This leads to five input fluxes, which may be thought of as J, the two net oxygen fluxes j l + j a at each surface, and 
the flux ratios j l /{j % + j a ). 

As a guide to solving the present problem, note that Ref. considered, in one dimension, slow steady discharge 
of a lead- acid cell for concentrations where reactions were unimportant. (Two and three dimensions have also been 
considered^) It was found that the mass loading (and for discharge, unloading), given by d t n (for ions H + and 
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HSO4 ), is uniform in space and time, so the background densities decreased uniformly in space and time. The 
continuity equation then implied that the fluxes vary linearly in space and are constant in time. It was then found 
that all of the vector quantities in this problem - two fluxes, two density gradients, and the electric field - vary linearly 
in space and are constant in time. The two ion densities have the same gradients, but they have a constant offset, 
leading to a constant non-zero charge density that is proportional to the current flow. This does not violate charge 
conservation because the screening modes at the surfaces can take up the excess charge. 

In the present case we assume that the bulk response has a reaction rate density r and ftn's that are uniform 
both in space and in time. As a consequence of JSJ, the charge density p also is uniform in space (d x p = 0) and 
in time (dtp = 0). For fixed J and j % + j a at each surface, but unspecified flux ratios j l + j a ), this leads to 
a consistent solution of the equations for this system; since the equations are linear, the solution is unique. The 
uniformity assumption leads to a local constraint on the fluxes, and thus specifies the flux ratios at the surfaces. By 
adding in appropriate amounts of the diffusion-reaction mode at each surface, the flux ratio of the total system can 
be made to correspond to any experimental values for the net flux ratios. 



1. General Considerations 



The charge continuity equation is given by 

d tP + d x J = 0. (36) 

The assumption that d t p — implies that J is constant in space and in time, consistent with our taking J to have 
the same value at each surface. Moreover, dtp = implies there is no electrical charging, so the electric charge and 
the field E and potential <p that it produces should all be constant in time. The assumption that d x p = implies, if 
p is non-zero, that the field varies linearly in space, and that the potential varies quadratically in space, results that 
we will show are consistent with the other equations. These properties are shared by the slowly discharging lead-acid 
cell and chemically loaded I^CuC^+a without reactionsi^ 

Adding the second two equations of Q yields what amounts to conservation of the sum of O 2- ions and O atoms: 

d t {n i +n a ) + d x {j i +j a ) = 0. (37) 

By <|37[1 . since the 9 t n's are assumed to be constant in space and time, the net oxygen flux jo — j l + j a varies linearly 
in space, and can be determined in the bulk from its values at the two ends of the sample. 

From the continuity equations (J7J), the constancy in space and time of r and the c^n's then implies that the slopes 
of all the fluxes are independent of space and time. Hence, in the absence of surface modes, these fluxes j a , j l and 
j are given by a linear interpolation between their values at the boundaries. We will see that, because the fluxes 
vary linearly in space, so do the other vector quantities in the problem: d x rii, d x n e , and E — E x — —d x (f). (We have 
already argued that E is linear in space.) 

Let us rewrite Poisson's equation (QJ in terms of E = —d x (f>. This gives Gauss's Law which, in one-dimension, reads 



d x E=-{-2n t + n h ). (38) 

e 



Taking the x-derivative of (|38|l . and using the flux equations to eliminate d x rii and d x nh, gives 



, , i dn h ^ , e,dm2j x dn h j h 



diE = -(4-li + -^)E + -(^^- - )• (39) 

e opt opt £ opt on dph ah 

This is solved by assuming that E is linear in x, so d^.E = 0. Then (|39|) gives E as a linear combination of the j's. 
Since the j's vary linearly in x, then so does E, and our assumption that E is linear in x is consistent. Explicitly, 
iJsSU gives 

E = ~M = g 1 g [-P± J - + 2p^] . (40) 
e(§^+4f^) L dp h a h dpi oti 

Since, unless the values of j' and j h at the surface conspire to give no linear term in E of (|40l) . E typically will vary 
linearly in space. Thus the potential <fi typically will vary quadratically in space. Moreover, d x E typically is uniform 
in space, so by l|38|l the charge density p — e(— 2n^ + rih) typically is uniform in space, consistent with our assumption 
that d x p = 0. 
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Application of d x p — to J5J gives 

2d x n t = d x n h . (41) 

Eas. ipUj) and l|4*T|) . and the relations iflU then give 

2d x rn = d x n h = - U a ^ 9 ^' )(2^-+4^-). (42) 

duu dm 

Since the j's are linear in space, so are these two density gradients. 

We have already assumed that r of l(T^|l is independent of x, or = d x r. Using (|4*T1) and l(T5|l leads to 



d x n a = -d x n % = dn * 8ni d x n t , (43) 

w a %±°, 

on a 



so all three density gradients are linear in space. 

As a consequence of (143() and 1421) . 1|19D yields the condition 

3 = ~a a ( 1 f-)d x n a = — j 4 + 2 — f. (44) 
on a ai ah 

Thus, for the bulk response, j a has a specific dependence on j 1 and j h . If (|44|l is not satisfied at either boundary, 
then diffusion-recombination modes are generated, with amplitudes determined in subsection IV. C. 

2. Specific Considerations 

We can now obtain explicit values for the three fluxes, the field, and the three density gradients. 

(a) With a knowledge of the mass flux jo = j 1 + j a at each end of the sample, of the current J = e(— 2j 4 + j' 1 ), 
and of the condition l|44|l . one can obtain explicit values for all three j's. This gives 

j 2 ^(J/e) 

f = — a " W ' (45) 

2 JO + (l + ^J/e)) 

J 1 + 2°, + ' V*"' 

ai a h 

(b) With a knowledge of the three j's and l|40|) one can obtain an explicit value for E. It is sufficiently complex 
and unilluminating that we do not present it. 

(c) With a knowledge of the three j's and E, by the flux equations H19|) we can obtain explicit values for the three 
density gradients d x n. Since the d x n's all are proportional to j a (as given above), we do not present them here. Like 
the three j's and E, the d x n's vary linearly in space. We measure the associated deviations Sn from x — 0, so that 



Sn = / d x ndx. (48) 



Let us now write each density n as the sum of its equilibrium value and its deviation An, where the latter 
consists of three terms: an offset term A that is constant both in space and in time (and has yet to be determined), 
a term at (where dtn = a) that is constant in space but linear in time (this corresponds to chemical loading at a 
uniform rate, and has yet to be determined; we must also satisfy at <C n^), and a spatially- varying term Sn that is 
constant in time but has a non-zero spatial variation (that has in principle been determined): 

n e = (n° h + A h ) +a h t + Sn h , m = (n? + A*) + atf + 6m, n a = (n° + A„) + a a t + 6n a . (49) 
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Note that nf — 2n°, since the system is neutral in equilibrium, and n° is (in principle) known from the equilibrium 
thermodynamics. 

We can now obtain explicit values for the three a's. 

(a) Since we assumed that dtp = 0, by we have 

2d t n l = d t n h . (50) 

Eq. JSUJ), when applied to i|l9*|l gives 

ah = 2ai. (51) 

(b) The requirement that r of Hlfijl satisfy dtr = 0, when combined with l|51J) . leads to 

a a = ah = jr, an- (52) 

w a 

on a 

(c) Since j % and j a are known by linearly interpolating their values at the boundaries, j t + j a is known. Hence (|37|) 
gives 

a h + a a = -d x (j i +j a ). (53) 

This, with (|51|l and l]52[l. then determines all three a's, thus specifying the time-dependences of the n's. 

We can now obtain explicit values for the three A's, which correspond to deviations from true equilibrium values. 

(a) Since E is known from (|4T])) . by so are d x E and p = ed x E. Thus, by JSJ, 

n h - 2n, = Ai - 2A h = t = S -d x E (54) 
e e 

is determined, where d x E is obtained from a version of H4U|) with the fluxes replaced by their (constant) gradients. 

(b) Eq. @ for dtnh can be used to determine r in terms of ah and d x j h . The three continuity equations have 
already been used twice, once for charge, and once for oxygen, so that this is their third and final allowable use. From 
r we can obtain A a by l|16|l . written as 

r « -(2whAh + WiAi - w a A a ). (55) 

(c) The arbitrariness in time-origin for the terms at in the n's permits us to set one of the A's to zero, so we may 
take 

A h = 0. (56) 

The relations l]54p. (|55|) . and (|56|l determine all three A's. 

Thus we have obtained the time variation at, the spatial variation Sn, and the offset A for each n. 

Note that the reaction rate density r, and the dn/dVs are independent of position and of time, as assumed. 

The steady-state bulk solution we have obtained is very general, having three free parameters (the current and the 
net mass fluxes at each surface), but the most general steady-state bulk solution has five free parameters. By adding 
in a diffusion-reaction mode at each surface, we can maintain the same current and net mass fluxes at each surface, 
yet still permit five free parameters. Hence, since the equations are linear, in obtaining one steady-state solution, we 
have obtained the steady-state solution. 

The results above hold even if the fluxes are slowly varying in time, provided one uses the instantaneous values of 
the fluxes. In this case the a's also are slowly varying with time, but their dominant dependence on time is given by 
the above approach. 



C. Using Surface Solutions to Match the Boundary Conditions 

Assuming arbitrary fluxes at the boundaries, we now determine how much of each type of response (bulk and 
surface) is generated. 

First, note that (I44|) applied to our specific bulk response (to which we append the subscript b), can be rewritten 

as 

jZ = ^f b + 2^j>\ (57) 
cti ah 
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Now note that the electric current density of (|24|) . 

J=e(j£-2jt), (58) 
is uniform in space, so it holds both in bulk (b) and at edges (E). With (|58[l . (|57|l can be rewritten as 

# =(4 *L+2V 6 + 2 --- ( 59 ) 
ai a/i e 

At either edge E, write 

= j| - Aj%, f b = j*. + Aj|, (60) 

where j% and j E are the (in principle known) total O and 2 ~ fluxes at the edge, and Aj% = — Aj E are the amplitudes 
due to the surface solutions. Substituting both parts of ij6TJ|) into then yields 

(1 + - + 4^)A?1 = 3% -(-+ *-)3h -2^ J -=j%- ^j E - 2^j E . (61) 
ah 014 au ai ah e a l ah 

Hence, from a knowledge of all three fluxes at an edge E, and the ratios a a /ah and a a /ah, the associated surface 
solution amplitude Aj E can be determined. By H29|) . Aj E /2 = Aj E = —Aj%. 

Global electroneutrality may be maintained by adding in screening modes, with the appropriate amplitude, at 
each surface. The specific value of the screening mode amplitudes will depend upon the nature of the contact at 
the surface, including the electrical contact resistance. Moreover, there can be a dipole layer at the surface, so the 
La2Cu04+5 itself need not be globally electroneutral, only the La2Cu04+,$ and a small region at each end of the 
adjacent materials. 

The results of this section are sufficiently general that they include those of Sect. IV. A, where there are no density 
gradients in the bulk, and uniform flow of ions (f) and holes (j h ). Note that if the surface values of f equal the bulk 
values, then no diffusion-reaction modes are generated at the surfaces. If the surface values of j i are zero (i.e. high 
pressure oxygen on one end and low-pressure oxygen on the other end), then the input flux of oxygen must come in 
the form of atoms with j a non-zero, and the diffusion-reaction modes then permit the system to convert atoms to 
ions within the diffusion-reaction length of the surface. 



V. ON COMPARISON TO EXPERIMENT 



Comparison to experiment requires a knowledge of the three thermodynamic derivatives d^i/dn, the three diffusiv- 
ities D (or the three a's), and the reaction constant A, which are independent of the fluxes applied to the system. 
Inclusion of the five independent fluxes at the surface means that the theory for steady state transport involves twelve 
independent quantities. The question then is how they may be determined. 

In principle, equilibrium measurements can yield the three dfi/dn's. Measurement of / for the screening mode can 
provide an additional constraint on the dfi/dri's. 

The fluxes and fields do not depend upon the reaction parameter A. It appears that A most easily can be obtained 
by measuring L for the diffusion-reaction mode, provided that the <9/^/(9n's already are known. 

Dc current measurements with electrodes that do not permit chemical loading, and ac current measurements with 
any set of electrodes, give cr, + av 

Measuring electrodes for holes and ions yield effective fields 

E h = — , Et = (62) 
ah 2a l 

Thus measurements with either type of blocking electrode (for holes or ions) gives two combinations of coefficients, 
one for jo = f + j a and one for J/e. Because the equations for the effective fields hold for all jo and J/e, these 
coefficients are overdetermined. 



VI. CONCLUSIONS 



We have developed the theory of slow steady transport for La2Cu04+a, including both loading and chemical 
reactions. We have employed the principles of irreversible thermodynamics, in which appear the thermodynamic 
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derivatives dfi/dn, the diffusion constants D = a(d^/dn), and the reaction constant A. There are two surface modes, 
one associated with screening, and the other associated both with diffusion and drift, and with chemical reactions. 
The screening mode is an equilibrium response, with no fluxes at all. The diffusion-reaction mode has no current flux 
or mass flux, but non-zero atom, ion, and hole flux. For a given uniform current density J and mass flux at each 
surface, the system has a special bulk response where the reaction rate density and the mass loading are uniform 
in space and in time. By adding in the diffusion-reaction modes at each surface, general steady-state solutions are 
obtained. 

As a consequence of the generality of the boundary conditions considered, the present work applies both to oxygen 
partial pressure loading^ and electrochemical loading? of I^CuCU+i. Moreover, given the complex nature of the bulk 
solution, it is highly unlikely that the material parameters are such that either oxygen partial pressure loading or 
electrochemical loading will be able to avoid generating diffusion-reaction modes. By varying the imposed current and 
mass fluxes at each surface, it should be possible to obtain a number of constraints on the parameters appearing in the 
theory. Thus it would be of great interest to be able to apply the present results to an actual system. Moreover, the 
diffusion-reaction may be responsible for the puzzling phenomenon of "turbulence of the voltage" near the electrodes, 
and measurement of its length L would further constrain the values of the material parameters. 
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